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1 ITERATIVE MMSE EQUALIZATION-DECODER SOFT 

2 INFORMATION EXCHANGE DECODING METHOD AND DEVICE 

3 FIELD OF THE INVENTION 

4 The present invention concerns the field of decoding of information symbols 

5 received over a channel, 

6 BACKGROUND OF THE INVENTION 

..^7 Information symbols are transmitted over channels in most information, storage 

^ and communication devices. The channel can comprise, for example, a wired or wireless 

^9 communication link between a transmitter or receiver, or an information path within a disk 

TO drive. More generally, any information path of limited bandwidth through which digital 

Jl signals are transmitted may be defined as a channel. Such channels are subject to various 

.^2 types of distortion, e.g., additive noise and intersymbol interference. These distortions limit 

^3 the data rate efficiencies for communication over such channels. Accordingly, much effort 

J4 has been devoted to reducing the effect of such distortions in the channels. The obvious end 

05 goal is to decode symbols accurately. More specifically, the goal is to have low decoding 

16 error rates. A second and conflicting goal is to increase the ratio of data symbols to total 

17 channel symbols transmitted, i.e., the data rate. 

1 8 Recent focus for joint equalization and decoding has been on iterative decoding 

19 schemes, and in particular, turbo coding schemes. These schemes are confidence building 

20 schemes in which an equalizer and decoder trade soft information in the form of symbol 

21 estimates until convergence is reached and hard decisions for symbols are output. These 

22 iterative schemes treat the channel contribution like an error correction code. The data 
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1 symbols are generally permuted before transmission and are typically protected with a 

2 convolutional error correcting code. On the decoding side, soft information is exchanged 

3 between an error correction decoder and a channel decoder. However, all current 

4 implementations of these soft information exchange schemes use a forward/backward, 

5 Viterbi, or similar decoding algorithm for the decoding fimction. Specific examples are the 

6 forward/backward and soft output Viterbi algorithm for convolutional codes. In such 

7 decoding schemes, the complexity of the decoder is a design choice, but the complexity of 

8 a trellis based equalizer is a fimction of the length of the channel impulse response. This 

9 limits practical usefiilness of the conventional decoding schemes to a class of channels 

10 having reasonably short impulse response lengths or small signal constellations. 

1 1 The complex nature of conventional equahzing schemes in soft input exchange 
Jj2 methods and devices also places certain physical limitations on devices which rely on such 
■^3 methods and devices for channel decoding. For example, there exist no solutions beyond one 
^ dimension for ft)rward/backward and Viterbi algorithms. This dictates a one dimensional 
Qp5 data stream, such as the stream obtained by a magnetic disk drive head or an optical disk 
46 drive. There is no significant physical reason for restriction to one-dimensional streams of 
37 data, but an as yet insurmountable hurdle for extension to higher dimensions is the lack of 
Jl8 higher dimensional forward/backward and Viterbi algorithms. A decoding method with 
39 multi-dimensional decoding capability would free devices fi*om the requirement of creating 

20 a one dimensional stream of data. For example, two dimensional regions of disk drives could 

21 be read simultaneously if two dimensional channel decoding is supported. The 

22 forward/backward and Viterbi algorithms used for data decoding are unable to handle the 

23 multidimensional task since multidimensional solutions are unknown. 

24 Accordingly, there is a need for an improved method for decoding data 

25 received over a potentially noisy channel which addresses some or all of the aforementioned 

26 drawbacks. It is an object of the invention to provide such an improved method. 
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1 SUMMARY OF THE INVENTION 

2 These and other needs and objects are met or exceeded by the present iterative 

3 MMSE equaUzation-decoder soft information exchange decoding method and device. The 

4 method uses a MMSE equaUzer which receives and outputs soft information. The equaHzer 

5 exchanges soft information with a soft information decoder. The simple nature of the 

6 equaHzer permits solutions beyond one-dimensional data streams and for channels of 

7 arbitrary length and for signal constellations of arbitrary size. 

8 BRIEF DESCRIPTION OF THE DRAWING 

9 Other features, objects and advantages of the invention will be apparent to 
10 artisans from the detailed description and the FIGURE, which is a block diagram of a data 
Jl communication system including an iterative MMSE equalization-decoder for soft 

information exchange decoding. 

m DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Jl 4 The invention is directed to a method for data recovery by confidence building. 

35 The invention is applicable, for example, to communications, telephony, recording, data 

56 transmission, and restoration. For example, even when the data has not been protected by 

9? an error correcting code, there may be other opportunities for confidence building decoders 

18 to operate. Suppose that the data were transmitted over multiple channels to the same 

19 receiver. This might arise when a preferred user in a cellular telephone system was given 

20 two or more cellular phone channels (for example multiple time-slots in a time-division 

21 multiple access scheme, multiple firequency slots in a frequency-division multiple access, or 

22 multiple spreading codes in a code division multiple access scheme), it might include storing 

23 the information on a data recording device multiple times (modular redundancy), or a 

24 receiver that has multiple antennas which can receive the same data at multiple locations. 
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1 According to the invention, even without error correction coding, a conjSdence building data 

2 recovery system may be constructed, in which equahzers for each of the different channels 

3 exchange their beliefs, or estimates, of the data, until a consensus is achieved. A typical and 

4 more general case arises when there is one or multiple channels over which data is sent, and 

5 the data is protected by an error correction code, likely different for each channel, before 

6 being sent over each channel. 

7 The generality of the confidence building data recovery method of the 

8 invention will be apparent to artisans. The method will be illustrated with reference to an 

9 error correction coded communication channel. An exemplary communication system 

10 including such a channel is shown in FIG. 1 . 

1 1 For simplicity, we consider transmission of binary digits b^ using a Binary 
& Phase Shift Keying (BPSK) modulation scheme. Artisans will understand that the invention 
^ is readily generalized to other modulation schemes. We assume that the bits bj are encoded 
M4 in a forward error correction code C (typically a convolutional code), by an encoder 10, and 
% that the bit stream is permuted by a suitable (e.g. random) interleaver 12 before the 
46 transmission. 

S? The encoded symbols are denoted aj. 11 is a permutation of the bits a ^ The 

Hb8 transmitted signal may be expressed as 

09 ^(0 = E b,h^it - iT)exp(jX2%t + (j)^)) 

i 

20 where 4 is the carrier frequency and 4>o is a constant angle. The received signal equals 

21 xit)=Y<b.h^it-iT)exp(jX2nf^t + ^^)) +w(0, 

22 where is a constant angle and w(t) is an additive noise term. In addition to intersymbol 

23 interference, a channel 16 typically also adds a noise component. The signal x(t) is 
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1 demodulated, passed through a receiving filter and sampled leading to the equivalent 

2 discrete-time model of the communication link which is depicted in FIG. 1 . 



3 x[u]=J^b.h[n-'i]^w[nl 

i 

4 where h[j] are the sampled values of the overall impulse response of the communication 

5 system. 

6 A discrete-time model of the channel 16 with intersymbol interference and 

7 additive noise Wj^ is described as 

% x[n\ =6„/2[0] +>v[«] b.h[n-i]. 

n9 where the term Y,i^n Ei hih[n - i] is usually referred to as intersjnnbol-interference. 
ift) We note that the received sampled signal x[n] is a noisy version the 

!l,l convolution of the bit-sequence with the overall impulse response. In particular, for a finite 

92 length impulse response h[j] the sequence x[n] can be represented as a noisy output of a 

ffS Markov process. The task of estimating symbols b^ in the absence of a coding constraint is 

34 referred to as equalization. The invention concerns joint decoding and equalization of the 

15 bit-stream x[n] in an iterative decoder 20 having a soft MMSE equalizer 22 exchanging 

16 information with a soft decoder 24. 

17 It is well known how to optimally estimate b^ from the sequence x[n] using the 

18 so-called forward/backward algorithm. However, the complexity of this approach is usually 

19 impractical because the complexity is essentially the product of the complexities of an 

20 optimal decoding procedure and an optimal equahzation only procedure. The complexity of 

21 an optimal equalization only procedure grows exponentially in the length of the overall 
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1 impulse response, which, for practical channels with impulse response lengths in the tens to 

2 hundreds, is infeasible= 

3 Iterative (Turbo) decoding and equalization uses at least two soft input/soft 

4 output (SISO) devices that communicate information. A SISO equalizer takes as input the 

5 received sequence x[n] and some additional information about the symbol sequence bj. This 

6 additional information could for example be a prior probability on the bits. The SISO 

7 equalizer produces a soft estimate of the symbols. This soft output could for example be a 

8 probability mass fimction indicating the likelihood of any symbol bj being equal to 0 or 1. 

9 A second SISO block performs the decoding of the estimated bit sequence with 

10 respect to the code. In other words the soft input decoder takes as input reliability 

1 1 information (typically a probability mass function) and produces as soft output a probability 
52 mass fimction indicating the likelihood of any symbol bj being equal to 0 or 1 . The input 
%3 reliability information for the equalizer (decoder) is denoted 'n:.^(7i.„). The output of the 

devices is denoted as ^out('^out) • overall channel impulse response h[j], let the 

^ received values x[n] and reliability vectors T^ij, [n] n = 1 ... N be given. Without loss of 

46 generality we assume that 0 < n^^ \j] < 1 for all j G {0 . . . N}. (The n^^ [j] can be thought of 

d|7 as prior probabilities). 

f|8 We define a SISO equalizer mathematically as a device that performs an 

S9 equalization fimction eq : tf^ x [0, 1]^ [0, 1]^ ^^(y^'^t)^^out^ where Tif^^ can be thought 

20 of as a vector of a posteriori probabilities. 

2 1 Similarly, we mathematically define an SISO decoder as a device that performs 

22 a decoding fimction dec: [0, 1]^ [0, 1]^ dec(7T;^)>^7i^^p where n^^^ can be thought of as 

23 a vector of a posteriori probabilities. 

24 The main idea of conventional iterative joint equalization and decoding is to 

25 let the two SISO blocks exchange information until they agree on a decision or until a 

26 maximal number of iterations is reached. As soon as the predetermined termination criterion 
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1 is satisfied, the output and the input of the SISO decoder are combined in order to form a soft 

2 combined equahzation decoding decision. If desired, a hard decision on any symbol can be 

3 derived from this soft value. 

4 Typically after a suitable number of iterations the output of the joint equalizer 

5 decoder is formed by multiplying the values of the likelihood ratios 71^/(1 -tt^) and 

7 In order to achieve good performance in such an iterative scheme, a 

8 permutation is included in the data path, which spreads out statistical dependencies between 

9 estimates. The exchange between decoder and equalizer establishes an inherent feedback 

10 loop. It is important that this feedback be mitigated as much as possible. We say that an 

11 equalizer (decoder) is well-behaved if tc^^^[/] is not a fimction of ti^-^O'] is not a 

12 fimction of 7t;^[/]). 

15 We now can describe the iterative equalization/decoding scheme in the 

^4 following algorithm: 

^ 1. INPUT: A permutation 11, a vector of received values y, and prior 

=\6 information about the symbols expressed as a vector of probabilities T^iJj^in usually a 

5? vector containing a value 0.5 in each position.) 

2. Repeat the following steps until a termination criterion is reached. 

§9 (a) T^lr = e?0,7T:f,) 

20 (b) < = n-k<j 

21 (c) C = deciTil) 

22 (d) < = n(C) 

23 (e) If the termination criterion is not satisfied go to step 2a. 

24 3 . OUTPUT hard decisions for symbol b^ comparing 

25 tol. (l-CWKlX™ 
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1 Iterative joint equalization and decoding schemes can be compared on the basis 

2 of the above algorithm and various combinations of decoding functions and equalization 

3 functions have been reported on. As decoding functions we explicitly mention the forward- 

4 backward and the SOVA algorithm for convolutional codes. Adaptations of these algorithms 

5 have been employed in the equalizer function. However, while the complexity of the decoder 

6 is a design issue, the complexity of a trellis based equalizer is based on the length of the 

7 channel impulse response. This severely hampers practical use of iterative equalization 

8 decoding techniques based on the turbo principle. For channels with long impulse responses, 

9 the exclusively used equalization techniques are either linear equalization or decision 

10 feedback techniques without use of the decoder in the equalization process. The invention 

1 1 solves this problem by providing an equalization algorithm independent of the length of the 
52 channel impulse response that can be used in iterative decoding. In particular, a minimxmi 

mean square error equalizer is demonstrated by the invention to facilitate joint decoding and 

a|4 equalization. The iterative technique of the invention is applicable to channels with a short 

55 or long impulse response, and with small or large signal constellations. 

46 Let b[n] be a sequence of symbols that have been transmitted through a 

fl|i7 communications channel, such as a wireless or wireline data link, or recorded onto some data 

fis recording medium, such as an optical ROM or magnetic disk or tape. A linear, equivalent 

29 base-band model for the received (or read-back) sequence of symbols is given by 

20 x{n\ = Kklbin -k] + w[w], 

21 where h[k] is the length + L2 + 1 impulse response of the channel, and w[n] is an additive 

22 noise term modeling electrical and thermal fluctuation and undmodeled components of the 

23 channel. This generic linear model is a widely accepted mathematical model for the 



-8- 



distortions induced on a sequence of symbols both in communications and data storage 
applications. A linear equalizer for this channel can be used to attempt to reduce the effects 
of the intersjmibol interference and the additive noise induced by the channel. A linear 
(affine) equalizer with coefficients c[n, k] and offset g[n] and symbol estimate b[n] can be 
expressed in the form 



where the equalizer coefficients are written as a function n to enable the possibility of 
different coefficients used to estimate each sjmibol ^ [n], and the offset g[n] provides a richer 
class of linear estimates which can account for a non-zero mean prior. The channel model 
can be written in matrix form as 



^2 

£ c[n,k]x[n+k]+g[n]. 



x[n] = HB[n] + w[n], 



where the channel response matrix, H, is given by 



h[-LJ 



0 



0 



0 



0 



0 



0 



h[L,] h[L,-\] ••• h[-L,] 



and the signal vectors are given by 

x[n] = [x[n - N, . . . x[n] . . . x[n + N^]]^ 
B[n] = [b[n - Ni - . . . b[n] . . . b[n + + 



1 w[n] = [w[n - NJ . . . w[n] . . . w[n + N2]]^. 

2 The equalizer output can also be expressed in matrix form, simply as 

3 6[n] = c[nrx[n] + g[n], 

4 where the equalizer coefficients c[n,k] are written 

5 B[n] = [c[n,-N,L...,c[n,N + 2]]. 

6 With this channel model, the mean-square error of a symbol estimate b[n] based on the 

7 estimate b[n] is given by 

8 E{\b[n]-b[n]\'}, (2) 

9 where the expectation is taken over the distribution of the symbols b[n] and the noise w[n]. 
10 In the conventional traditional approach to the design of complexity-constrained (finite Ni 
M and N2) minimum mean-square error (MMSE) linear equalizers, it is assumed that the 
jp symbols b[n] are equally likely to take on all possible symbol values, and that there is no 
r© additional information about their values available. The conventional equalizer is determined 
^4 by finding the coefficient values c[n] and g[n] which minimize the mean squared error (2), 
^jL5 which, since the symbols are assumed unknown and equally likely for all time n leads to a 
ffl6 single set of coefficients, c, and the offset is given by g[n] = 0. 

rlf7 According to the invention, the equalizer 22 implements an equalization 

58 algorithm which has available a set of priors over the symbols. For example, if the symbol 

19 alphabet is binary, then this would correspond to the availability of the sequence ^iJ.n] = 

20 Prob {b[n]==l } . In the sequel, we assume that the channel response h[n] is real and that the 

21 symbol alphabet is b[n] e {-1, 1} for simplicity. Extension to complex baseband channels 

22 and higher-order symbol constellations is straightforward. In this case, the MMSE equalizer 

23 22 can be designed incorporating these priors into the optimization. Hence, the equalizer 

24 coefficients, c[n] and g[n] can be determined by finding the minimum of the mean-squared 

25 error (2), where the expectation in (2) is over both the additive noise in the channel, and the 
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given (time-varying) prior over the symbols. As a result, the equalizer coefficients will vary 
with time index, n. This leads to the following formulation, 
b[n] = E{b[n]}+[E{b[n]B[n]}H'^ - E{b[n]}E{B[n]^}HT 

[HE{B[n]B[nf }H^ + E{ | w[n] P}I - HE{B[n]}E{B[nf }H^-^ (x[n] - HE{B[n]}). 
Once the equalizer has produced MMSE linear estimates of the symbols b[n], these 
estimates must be mapped into priors n^^^ . One method for mapping the outputs of the linear 
equalizers is to assume the output distribution b[n] is conditionally Gaussian, distributed 
about the symbol values. This leads to the following mapping 



Prob{b[n]=l\b[n]} 

2 



1 -^tanh 



( \ 
b[n] 



\ 



\ 



where a. is the variance of the conditional output distribution given the symbol b{n\ = sign (6 [w]) . 
In order for this equalizer to be considered well-behaved, the estimate 6 [n] cannot be a 
function of [n]. Hence, the expectations must be taken over a distribution of the symbols 
which excludes 7T:^[n] for the calculation of 6[n]. However, in calculating 6 [k], k n, 
7r-„[n] may be used. This leads to the following method for computing the output 

E 

distribution given the observations, x[n] and the input distribution iij^. 



CREATE BUFFERS 

1 . Create buffers for the priors, the signal x[n], the expectations M 
E{b[n]B[n]}, the correlation matrix B[n] = E{B[n]B[n]^}, and the means [n] = E{B[n]} 

7t(n) A [7t:(")[-Ni - L2], - . . , 7i:(">[0], . . ., 7t:(")[N2 + Li]]'^ 

jt(n) A [X(")[-Ni], . . X(")[0], . . ., XW[N2]]^ 

bb A[bb(")[-N, - LJ, . . ., bb(")[0], . . . bb<")[N2]]T = [0, . . . 0, 1, 0, . . 0]T 



1 2. Initialize buffers for priors fc^") and data x^"\ in terms of the signal x[n] 

2 and the input ti^. 

3 x(«> = [0,0,,..,x[0],x[lL...,x[N2]r 

4 7t(«) - [0, 0, . . 0, n%[Ol n%[ll . . . , ^[N^ + L,]f 

5 3, Loop over the data for n== 0, . . N: 

6 71^^) [0] - Vi 

7 mb^''^ = 2ii^^^-l 

^ ^ mb mo 

9 diag(B) = diag(l,l,...,l) 

10 c[n] = [H (B - ^-^6 ^"^mb ^"^'0 + a ^I] H ^l, 

U 6 [n] = ^-fi + c(")^ (x(") - H ^-^z, (")) 

m x("^'^ = [xW [-Ni + 1], . . . , x("> [N2], 0] 

P fc("+i) = [tt(") [-Ni - L2 + 1], . . . , 7t("> [N2 + L,], 0] 

p ifn<N-N2 

?5 x("+'>[N2] = x[n + 1 + N2] 

=16 ifn<N-N2-L, 

pP7 7i:("+i>[N2 + Li] = Kf^ [n + 1 + N2 + Li] 

318 4. Estimate output variance al = (yar(b \b>0)+ var(b | b<0))/2 

r|;9 5, Determine output priors, Ti^^y, = 1/2(1+ ^awA (— ^)) 

20 SOFT-INPUT SOFT-OUTPUT DECISION FEEDBACK EQUALIZATION 

21 For a minimum mean-square error (MMSE) decision feedback equalizer, the 

22 channel model can be written in similar matrix form to the MMSE linear equalizer, 
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where, for simplicity of notation, it is assumed that the number of decision symbols feedback 
to the equalizer is given by M = Ni + L2, for 3[n] = [d[n-l], . . d[n -M]]. For binary 
antipodal signaling, d[n-k] = sign (6 [n-k]), and for higher order signaling constellations, a 
suitable quantizer to the nearest symbol would be used. 

The MMSE DFE coefficients are then given by 

b[n] =E{bln]} + 



E{b[n]b[n]}H^ -E{b[n]E{b[nf}H''\ 0 



UM 



HE{Un]blr'f)H^^E{w[nf}I-HECb[n]}E{i[nY}H ^g£{g[«M«f } -HE{lln]}EQ[nf } 



3c W -HE{b[n]} 
dW -E{d[n]} 



When this MMSE DFE is made well-behaved, we set E{b[n]} = 0, which, together with 
some algebra, reduces this expression considerably, to 

b[n] =E{blnMfif}H ^ [HE{i[nMnf}H ^ +E{-w[nf}I -HE{b[n\}ECb{nY}H ^ - 
(HE{mdlnf}-HE{l[n]}E{2nf})E{2lnMn]T\EQln]^^^^ 
[C[n] -HE{b[n]) -(HE{b[nMnf} -HE{m}E{d[rif})EQ[nMnfy\dW -EQin]}) . 
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1 The preferred implementation of the present invention is shown in the 

2 FIGURE. In the FIGURE, a process for encoding/transmission and reception/decoding 

3 begins with a set of digital data, depicted as data bits aj. These data bits are then encoded, 

4 using forward error correction coding 10, to produce the encoded sequence of bits, depicted 

5 as a.. 

6 The encoded data bits are now interleaved (re-ordered) in time using the data 

7 interleaver 12. The purpose of the data interleaver is to re-order the data such that the 

8 statistical dependencies between the data bits are spread out in time. This makes adjacent 

9 data bits in the re-ordered sequence, bj in the figure, approximately independent of one 

10 another. The re-ordered bits bj are now ready to be transmitted over the channel 16. The 

1 1 process of mapping the bits b^ into channel symbols and transmitting them over the channel 
is depicted in the FIGURE as block 16. The ISI channel 16 introduces distortion into the 

J3 sequence of channel symbols. The channel 16 also is assumed to exhibit additive noise, as 

J|4 depicted in the FIGURE. The output of the channel 16, is then the sequence of corrupted 

'"f5 channel symbols, x[n]. The received sequence x[n] is then processed by the receiver block 

!l^6 20 to remove effects of the channel. The receiver block comprises several elements. 

W First, the received data symbols are equalized using a soft-input/soft-output 

fflS MMSE equalizer 22. The equalizer 22 attempts to eliminate the intersymbol interference 

59 (ISI) from the channel. The output of the equalizer 22 is a set of priors, or confidence levels, 

20 in the symbol values, labeled as II in the figure. The two inputs to the SISO MMSE 

21 equalizer 22 are the channel output symbols x[n] and a set of confidence levels in their 

22 values, II , On the first pass through the equalizer, an initialization set of confidence levels 

23 are used, which are equally-likely to take on all values. Subsequent passes use confidence 

24 levels produced by the SISO decoder 24. 

25 The confidence levels are then re-ordered using the de-interleaver 12a to place 

26 them in the same order as the corresponding bits in the encoded sequence a^. The confidence 
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1 levels can now be used in a soft-input/soft-output decoder 24 to produce estimates of the 

2 original uncoded data bits a,, labeled as II out figure. The confidence levels over the 

3 sequence aj are then interleaved 12 again, back to the ordering of the channel symbols and 

4 the interleaved data bits bj. The confidence levels are now used as input to the SISO MMSE 

5 equalizer 22, together with the data x[n]. Whereas in the first pass through the equalizer 22, 

6 the confidence levels were arbitrarily preset to initialized values, now the confidence levels 

7 have been determined by the decoding process. This cycle is repeated until either a 

8 convergence criteria is met, or a sufficient number of passes over the data have elapsed. 

9 Typically, the convergence criteria will consist of a prespecified measure of match between 

10 the confidence levels determined by the decoding process and those determined by the 

1 1 equalization process. Other possible convergence criteria could include testing that the 
M confidence levels determined by one or the two SISO devices have not changed appreciably 
J3 over a sequence of passes. 

3fl While various embodiments of the present invention have been shown and 

%^ described, it should be understood that other modifications, substitutions and altematives are 

46 apparent to one of ordinary skill in the art. Such modifications, substitutions and altematives 

£^ can be made without departing firom the spirit and scope of the invention, which should be 

di8 determined fi^om the appended claims. 

f|i9 Various features of the invention are set forth in the appended claims. 
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WHAT IS CLAIMED IS: 



1 LA method for estimating digital data received over a potentially noisy 

2 channel which adds intersymbol interference or additive noise, or a combination of 

3 intersymbol interference or additive noise, the method comprising steps of: 

4 inputting data received from the noisy channel into a SISO MMSE equalizer; 

5 inputting a set of priors over symbol values of the noisy channel including a 

6 separate prior for each received noisy channel symbol value, into the SISO MMSE equalizer; 

7 equalizing, by an MMSE equahzation in the SISO MMSE equalizer, the data 

8 received from the noisy channel and the set of priors over s)aTibol values to produce a symbol 

9 value estimate; 

H) mapping output of the SISO MMSE equalizer onto priors over the symbol 

%i values to produce a confidence indication in each of the symbol value estimates as a fimction 

22 of time. 

^ 1 2. The method according to claim 1, ftirther comprising a step of setting 

Q}2 parameters of the SISO MMSE equalizer according to MMSE crieterion over statistics of the 

ff3 channel noise and statistics of the symbol values. 

1 3 . The method according to claim 1 , wherein digital data transmitted over 

2 said potentially noisy channel is error correction encoded prior to transmission, the step of 

3 mapping output of the SISO MMSE equalizer comprising steps of: 

4 passing output of the SISO MMSE equalizer into a SISO error correction 

5 decoder; 

6 using an output of said SISO error correction decoder as the set of priors over 

7 symbol values; 
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repeating all steps of the method until a predetermined convergence criterion 
is reached between said SISO error correction decoder and said SISO MMSE equalizer. 

4. The method according to claim 3, wherein digital data transmitted over 
said potentially noisy channel is interleaved prior to transmission and a step of de- 
interleaving is conducted on the output of the SISO MMSE equalizer prior to said step of 
passing output of the SISO error correction decoder is interleaved prior to said step of 
repeating. 

5. The method acording to claim 4, wherein said step of equalizing 
excludes symbol value estimates which are functions of an input distribution of a current 
symbol being equalized. 

6. The method according to claim 5, wherein said SISO error correction 
decoder has its output restricted to exclude symbol value estimates which are functions of 
an input distribution of a current symbol being decoded. 

7. The method according to claim 1, wherein said step of equalizing 
comprises a fast update equalization of the order M^ (quadratic in the number of parameters 
M) which exploits redundant computations in successive equalizer computations. 

8. The method according to claim 7, wherein the fast update equalization 
is performed by applying the matrix inversion lemma to a matrix to be inverted in a design 
of equalization coefficients for the SISO MMSE equalizer. 
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1 9. The method according to claim 1, wherein said data is multi- 

2 dimensional and is error correction encoded then interleaved prior to transmission in the 

3 channel, and wherein said SISO MMSE equalizer handles single-input multiple output data, 

4 the mapping step comprising steps of: 

5 de-interleaving outputs of the SISO MMSE equalizer and re-serializing data 

6 output from the SISO MMSE equalizer into one-dimensional encoding; 

7 SISO decoding the re-serialized output of the SISO MMSE equalizer; and 

8 repeating said equalizing and mapping steps until a predetermined convergence 

9 criterion is met between said equalizing and SISO decoding steps. 

11 10. A method for equalization and decoding of digital data received over 

M2 a multiple noisy channels which each add intersymbol interference or additive noise, or a 

Q3 combination of intersymbol interference or additive noise, wherein data transmitted over each 

BE 4 channel is interleaved prior to transmission, the method comprising steps of: 
^5 performing a soft equalization for each channel by 

^^6 inputting data received from said one of the noisy channels into a SISO MMSE 

equalizer; 

8 inputting a set of priors over symbol values of the noisy channel including a 

9 separate prior for each received noisy channel symbol value, into the SISO MMSE equalizer; 

10 equalizing, by an MMSE equalization in the SISO MMSE equalizer, the data 

1 1 received from the noisy channel and the set of priors over symbol values; and 

12 mapping output of the SISO MMSE equalizer onto priors over the symbol 

13 values to produce a confidence indication in each of the symbol values as a ftinction of time, 

14 then, using de-interleaved output from said mapping step for one channel, iteratively 

1 5 decoding information for a second channel by repeating said performing a soft equalization 
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while substituting said de-interleaved output for said set of priors and substituting output 
from said mapping step for said second channel for said data received until a predetermined 
convergence criterion is reached. 



11. A data decoding device comprising: 
a SISO MMSE equalizer; 

a SISO decoder, the decoder exchanging symbol estimates with the SISO 
MMSE equalizer, the SISO MMSE equalizer produces MMSE linear estimates of transmitted 
symbols bXn] , and mapping the linear estimates to an output set of priors over the symbols 

E 



12. The device according to claim 11, wherein said equalizer maps the 
estimates by treating the output distribution b[n] as conditionally Gaussian, and distributed 
about the symbol values. 



13. The device according to claim 12, wherein the output distribution 
mapping is defined as: 



Prob{b[n]=l\b[n]} =- 

2 



1 + tank 



{ \ 
h\n\ 



where is the variance of the conditional output distribution given the symbol b\n\ = sign (6 [w]). 
and the estimate 6 [n] cannot be a function of 7i;^[n], and expectations are taken over a 
distribution of the symbols which excludes 7T;^[n] for the calculation of 6 [n]. 
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1 14. The device according to claim 12, wherein the following steps for 

2 computing the output distribution given the observations, x[n] and the input distribution Ujj^ 

3 is used by the MMSE equalizer: 

4 a. Create buffers for the priors, the signal x[n], the expectations bb M = 

5 E{b[n]B[n]}, the correlation matrix B[n] = E{B[n]B[n]^}, and the means ^b M ^ E{B[n]} 

6 71^") A [7i(")[-Ni - L2I . . . , 7i(")[0], . . 7C(">[N2 + Li]]^ 

7 x^"^ A [x^">[-Ni], . . x(">[0], . . , x^")[N2]]^ 

8 A[bb^")[-N, - L2I . . bb(")[0], . . . bb(")[NJ]^ [0, . . . 0, 1, 0, . . 0]^ 

9 b. Initialize buffers for priors n and data x^"\ in terms of the signal x[n] 
10 and the input tt^. 

n x('^ = [0,0,...,x[0],x[l],...,x[N2F 

iz 7i(^> = [0, 0, . . 0, 7r^[0], t4[1], . . . , <[N2 + U]Y 

c. Loop over the data for n= 0, . . ., N: 

Ik n^"'^ [0] = Vi 

|5 ;^(") = 27t^>-l 

|7 diag(B) = diag(l,l,...,l) 

|8 m = [H (B - ^<"U^"^') + o ^I] H 

m 6 [n] = ^("> + c(")^ (x(") - H ^(")) 

20 x("+*> = [x(") [-Ni + 1], . . . , x(") [N2], 0] 

21 Tt("^') = [7i("> [-Ni - L2 + 1], . . . , 71^") [N2 + Li], 0] 

22 ifn<N-N2 

23 x('^i>[N2]=x[n+l+N2] 

24 ifn<N-N2-Li 

25 7x("+i)[N2 + L,] = Ttfj^ [n + 1 + N2 + L,] 

26 d. Estimate output variance o^g = (yarib\b>0) +var(6 16<0))/2 
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27 e. Determine output priors, n^jjj, = 1/2(1+ tank (-^^)) . 
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1 ITERATIVE MMSE EQUALIZATION-DECODER SOFT 

2 INFORMATION EXCHANGE DECODING METHOD AND DEVICE 

3 ABSTRACT OF THE DISCLOSURE 

4 A method and device employing an iterative MMSE equalization-decoder soft 

5 information exchange decoding. The method uses a MMSE equalizer which receives and 

6 outputs soft information. The equalizer exchanges soft information with a soft input soft 

7 ou^ut decoder, preferably an error correction decoder. The nature of the equalizer permits 

8 solutions beyond one-dimensional data streams and permits solutions for long channel 

9 lengths and multi-dimensional data since the solution is not a function of the channel impulse 
10 response. 

Zl F \DATA\WP60\1201\63069\63069FIN APP 
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